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1 Introduction 


When flight first began, pilots controlled aircraft through direct force. The pilots moved control 
sticks and rudder pedals linked to cables and pushrods that physically moved the control surfaces such as 
the wings and tail of the plane. However, as power and speed of flight increased, more force was needed 
to control the aircraft and therefore hydraulically boosted controls were incorporated into the aircraft. In 
the 1960's, the idea of flying aircraft with electronic flight control systems was introduced. Wires 
replaced cables and pushrods, which in turn gave the aircraft designers greater flexibility in the size and 
placement of components to control the aircraft. A fly-by-wire system would be smaller, more reliable, 
and in military aircraft the systems would be much less vulnerable to battle damage [1], A fly-by- wire 
aircraft would be more responsive to pilot control inputs with the results being improved performance and 
design of a more efficient, safer aircraft. The quality of flight was greatly improved with this 
instantaneous sensing of pilot inputs. 

Fly-by-wire systems are safer because of their redundancies. They are more maneuverable 
because computers can command more frequent adjustments than a human can. Fly-by-wire is also more 
cost efficient because it is lighter and takes up less space than the hydraulic system, which in turn either 
reduces the amount of fuel needed or increases the amount of passengers/cargo the aircraft can carry [2], 

However, with this new concept of fly-by-wire comes awareness that the system, which now uses 
individual wires instead of cable bundles, is also more vulnerable to certain environmental conditions. 
Planes flying through adverse operating environments experience a phenomenon known as 
electromagnetic interference (EMI). There are many factors, man-made and natural, that can contribute 
to these phenomena including [3]: 
radar, 

- lightning, 

AM/FM/TV broadcast stations, 

industrial, scientific, and medical (ISM) equipment, 

automobile ignitions, 

personnel electrostatic discharge, 

esoteric nuclear electromagnetic pulse (NEMP), and 

power supply noise and switching transients inside electronic equipment. 

The electromagnetic interference phenomena dealing with radio frequency (RF) wavelengths, 
namely radar and AM/FM/TV broadcast stations, have become known as HIRF, High Intensity Radiated 
Fields. HIRF is a non-ionizing electromagnetic energy that is external to the aircraft. HIRF can cause 
adverse effects to the electronic equipment onboard the aircraft, which in turn may affect the safety of 
flight and landing [4], Electromagnetic fields may cause electrical signals to be induced on the aircraft's 
wiring and these signals can propagate to other electronic equipment that may cause a functional error 
known as upset [5]. 

Upset phenomena that can be caused by electromagnetically induced signals include [6] - [11]: 
change in data values of the input/output circuitry, 

logic changes on the data bus, address bus, and control lines of the processors, 
logic changes in registers of the central processing unit (CPU) of the processors, and 
logic changes in the arithmetic logic unit (AFU) within the CPU of the processors. 
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Upset phenomena such as these can interfere with normal operation of the processors within a control 
computer and result in control law calculation errors that can affect performance and reliability at the 
closed-loop system level [5]. 

Aircraft systems, critical and essential, are vulnerable to atmospheric electricity hazards [12]. 
The number of electrical/electronic systems aboard an aircraft is increasing. These systems are 
vulnerable to electromagnetic fields caused by HIRF and need to be tested to ensure they are not 
susceptible to HIRF. Another potential problem can occur if the aircraft is struck by lightning. The 
associated electromagnetic field may cause voltage and/or current transients to be induced into the 
electronic equipment. These transients can be produced in two different ways: the aircraft's interior may 
be penetrated by the electromagnetic field or the structural voltage may rise due to current flow on the 
aircraft. These are referred to as indirect effects because they may not physically damage the aircraft 
[ 12 ]. 


Electromagnetic fields can also penetrate inside the electronic equipment through imperfect 
seams, leaky connector apertures, and cracks in the protective shielding. Due to the increase of reliance 
on the electronic equipment of an aircraft for flight, adequate susceptibility methods must be utilized to 
ensure safe flight. Another factor that needs to be considered for ensuring flight safety is the possibility 
of reduced electromagnetic shielding by replacing the aircraft's metal skin with one made of composite 
materials. 

The goal of this research is to develop and demonstrate technologies that would eliminate vehicle 
system and/or component malfunctions as a factor in aviation accidents. In order to move closer to this 
goal, a mathematical model of electromagnetic fields coupling into the test equipment will be developed 
and validated. This model will then be used to extrapolate/predict measurements outside of the data set in 
order to determine the level of field strength that can be applied without damaging the test equipment. In 
Section 2 the open-loop calibration experiments and the resulting data are described. In Section 3 a 
statistical analysis of the open-loop data is explained. Section 4 contains the conclusion of the research. 


2 Description of Open-Loop Experiments 

The Systems and Airframe Failure Emulation Testing and Integration (SAFETI) Laboratory at the 
National Aeronautics and Space Administration (NASA) Langley Research Center (LaRC) is used to 
study the effects of HIRF on complex avionic systems and control system components. Linked with the 
High Intensity Radiated Fields Laboratory, also at LaRC, tests are being conducted on a quad-redundant 
fault tolerant flight control computer to establish upset characteristics of an avionics system in an 
electromagnetic field. This section of the paper describes the open-loop calibration experiments and the 
data collected [13]. 

The block diagram in Figure 1 represents the equipment used to collect the open-loop experiment 
data [14]. The flight control computer consists of four independent processors. Each processor has a 
1750 processor, 48 Kbytes of EPROM (Erasable Programmable Read Only Memory), 2 Kbytes of 
scratchpad RAM (Random Access memory), and 8 Kbytes of sharable RAM. All input/output is memory 
mapped into a 2 Kbyte space. Although digital in design, the flight control computer receives sensor 
inputs and sends actuator command outputs via analog voltages. These signals are interfaced to the flight 
control computer via nine shielded cable bundles that are approximately eleven feet in length. 


2 



Reverberation Chamber 



Figure 1 . Equipment used to collect the open-loop data. 


In addition, each line has a passive filter at the connection point to reduce noise on the signal 
lines. All nine cable bundles are attached to the analog electro/optic converter, which is a custom-built 
circuit that converts all of the analog/discrete inputs and outputs from the flight control computer into 
light that can be transmitted down nine single -fiber optic lines. This permits a safe, noise-free method of 
transmitting the signals in and out of the electromagnetic test chamber. 

In addition to the analog lines, a MIL-STD-1553 interface can be used. These four coaxial lines 
are bussed into a custom-built fiber optic converter. This interface can be used for digital input/output to 
close the control loop. 

A picture of the equipment set up for the open-loop experiments is shown in Figure 2. The flight 
control computer, also known as the equipment under test (EUT), is placed on a styrofoam block in the 
middle of the chamber. The large paddle on the left stirs the electromagnetic fields to obtain a 
statistically near homogeneous radiation environment during testing. The probes that measure the field 
strength inside the chamber are on stands to the left and right of the flight control computer. The source 
of power and signals are transferred through the cables that are connected to the bulkhead, which is the 
panel on the right-hand wall. The antenna on the right emits the radiation into the chamber through one 
of the cables connected to the bulkhead. The equipment that controls the radiation environment is located 
in a separate control room outside of the test chamber. The signals are passed to the larger box on the 
floor, which is the analog electro/optic converter. This is where the signals are converted to optical 
signals to be passed through the bulkhead. The smaller box on top of the converter is the 1553 box. 
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The flight simulator hardware is based on a twenty-slot VME backplane, Figure 3. The system 
consists of three Motorola 68040 based real time controllers, five digital to analog converters (DAC), one 
analog to digital converter (ADC), and one MIL-STD-1553 interface board. The five DAC’s and one 
ADC link the simulator to the analog electro/optic converter, and the MIL-STD-1553 board links to the 
digital electro/optic converter. Even with five DAC’s and one ADC, only one set of input/output lines 
can be supported. Consequently, the analog electro/optic converter performs a fan out function so that all 
four of the flight control computer’s processors receive the same inputs. 

The electromagnetic test chamber shown in Figure 2 is a 13 x 23 x 9 x /i foot mode-stirred 
reverberation chamber located in the High Intensity Radiated Fields Laboratory at LaRC [15]. Figure 4 
shows a picture of the outside of the test chambers. These are enclosed steel rooms that have been 
validated for shielding effectiveness to 120 dB to ensure that no radiation leaks outside of the chamber. 
The door is pneumatic in nature and has a seal that expands after the door is shut to ensure proper 
shielding. In essence, this is like a large microwave oven that provides a near homogeneous radiation 
environment. Using this type of test chamber enables the flight control computer to be exposed equally 
from all angles, so that the angle of incidence for maximum susceptibility does not need to be found or 
assumed. 



Figure 2. Flight control computer placed in reverberation chamber to collect data. 
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Figure 3. Flight simulation hardware. 



Figure 4. Outside look at reverberation chambers. 
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For the open-loop experiments, two special conditions were implemented. First, the flight 
simulation software was removed and the flight simulator was reprogrammed with a special calibration 
procedure. This procedure sent pre-selected reference voltages to each processor in the flight control 
computer in order to measure the noise level on the analog lines. The received voltages were then stored 
for comparison to the sent reference voltage. 

The second condition was that, with the above calibration procedure in place, the passive filters 
were selectively removed to introduce electromagnetic energy in the cables. This allowed the long 
electrical signal lines to couple electromagnetic energy into the flight control computer. These filters are 
extremely efficient at keeping the flight control computer radiation tight, so by removing specific 
combinations of filters the flight control computer processor could be weakened to the radiated energy. 
These combinations are shown in Table 1. The ninth cable cannot be removed because this supplies 
power to the flight control computer. 


Table 1. Passive Filters Selectively Removed for the Open-Loop Experiments 


Filters Removed 

Test Objective 

1,5 

Weaken Processor 1 

2,6 

Weaken Processor 2 

3,7 

Weaken Processor 3 

4, 8 

Weaken Processor 4 

1,2, 3,4 

Upper connector matrix cross-talk effects 

5, 6, 7, 8 

Lower connector matrix cross-talk effects 


The maximum and minimum voltages were calculated for each analog signal line and three 
critical points were defined: fifty percent of the maximum voltage, the zero point, and fifty percent of the 
minimum voltage, as shown in Table 2. Appendix A lists the actual reference voltages sent for each 
signal. 


Table 2. Pre-Selected Reference Voltages Sent to the Flight Control Computer 


Voltage Level 

+15 to -15 v 

+10 to 0 volts 

50% of Max 

+7.5 

+5 

Zero 

0 

0 

50% of Min 

-7.5 

0 


One hundred percent of the maximum and minimum values were never sent to the controller 
because with the added energy from the electromagnetic fields, the cumulative effect could have 
exceeded the hardware’s specification and resulted in component damage. Table 3 shows the effect the 
passive filters have on the system. The values shown reflect the overall maximum and minimum voltages 
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collected during testing. When the passive filters remain on the cables, there may be a nominal amount of 
radiation that enters the system that could minimally affect the voltage values received, see Table 3a. 
However, when the passive filters are removed, there is a change in the data collected, see Table 3b, 
which could have led to equipment damage if the maximum voltage threshold had been exceeded. 


Table 3(a). Voltage Values Collected with 
Passive Filter On 


Sent 

Received 



Min 

Max 

-4.995 

-5.034 

-4.832 

0.000 

-0.063 

0.182 

4.990 

4.958 

5.153 


Table 3(b). Voltage Values Collected with 


Passive Filter Off 


Sent 

Received 



Min 

Max 

-4.995 

-5.230 

-0.859 

0.000 

-0.216 

4.113 

4.990 

4.790 

9.035 


Figure 5 shows the flowchart for the automated calibration procedure. For a single run of the 
program 525,000 voltages were collected. The program was run for each of the six filter combinations 
listed in Table 1 which nets a total of 3,150,000 data points per test. 



Yes 

▼ 


Store 3,000 samples x 
35 lines x 5 voltage 
values as a file. 



Figure 5. Flowchart for specialized calibration procedure. 


Once the nominal characterization data was collected, a test plan was developed to obtain data 
from the flight control computer on the effects of radiated fields within the test chamber. Testing started 
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at 250 MHz and was incremented by 100 MHz steps until reaching a maximum of 950 MHz. Field 
strengths were incremented as shown in Table 4. 

The procedure followed for testing was to start at the lowest frequency and continue increasing 
the field strength to the highest level in Table 4. each time collecting the 3,150,000 data points described 
previously. However, tests were not conducted at all frequency/field strength pairs. Entries in the test 
matrix of Figure 4 indicate the test was halted if a large number of observed voltages violated a 
predetermined threshold, typically one percent above normal system drift. This test procedure was used 
for the indicated test configurations referenced in Table 1. 


Table 4. Open-Loop Experiment Test Matrix 


250 

50 

100 

150 

Fie] 

200 

Id Strer 

300 

igth (v/i 

400 

n) 

450 

500 

550 

600 

X 

X 

X 

X 







350 

X 

X 

X 

X 







450 


X 


X 

X 

X 

X 

X 

X 

X 

Frequency 55 q 


X 


X 

X 

X 

X 

X 

X 

X 

( MHz ) 650 


X 


X 

X 

X 

X 

X 

X 

X 

750 


X 


X 

X 

X 

X 

X 

X 

X 

850 


X 


X 

X 

X 

X 

X 

X 

X 

950 


X 


X 

X 

X 

X 

X 

X 

X 


Data were stored in an ASCII file with five columns, one for the reference voltage sent and one 
column for voltages collected at each of the four processors. Each data file included three thousand 
samples, one thousand samples for each of the three reference voltage groups. Figure 6 shows a partial 
data file collected. Due to the overwhelming amount of data collected, one of the first tasks was to 
develop a data management scheme. This scheme involved standardizing file conventions such as file 
naming, data format, and data variables collected from the tests. It also involved storing the data to a 
medium that was accessible to all users and creating web pages that listed the unique test parameters and 
the directory locations of all files. 


3 Statistical Modeling of Open-Loop Experiments 

In order to examine the open-loop test data, a reduced representation was required. The first 
effort involved plotting the data for all frequency/field strength combinations. Figure 7 shows an 
arbitrarily chosen plot of the test data for “spl”, left spoiler in degree. Each subplot represents a different 
set of data for the signal. The columns depict the reference voltage sent while the rows are the field 
strengths. The first column of subplots represents fifty percent of the minimum voltage sent, the second 
column represents the zero point, and the third column represents fifty percent of the maximum voltage. 
For comparison sake, all column axes have the same values. 

The first row shows the three reference voltage levels at one field strength. Each subsequent row 
shows an increasing field strength at the same reference voltage level as that above it. As the field 
strength increases, the plots flatten out and the data begins to spread over a broader voltage range. This 
trend appeared in all of the test data examined. 































































Next, distribution curves were fit to the data since the data had no well-defined shape. The 
distributions that were fit to the data include the Uniform, Normal, Weibull, Rayleigh, Beta, Gamma, 
Exponential, and Pareto distributions [16, 17], The distributions fit to the data, along with their 
associated probability density functions, are shown in Table 5. Using the method of least squares, it was 
determined that the curve that best fit the data overall was produced by the Normal distribution and 
therefore from this point forward this work will concentrate on the Normal distribution. The equation 
used to define the Normal distribution, as shown in Table 5, has x as the voltage value, p the mean of the 
data, and a 2 the variance of the data. 
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Figure 6. Partial open-loop experiment data file collected. 

Figure 8 shows the same data as Figure 7 with the Normal distribution curve overlaying the data 
values. For this particular signal, the Normal distribution has a good fit. There were some signals in 
which the Normal distribution did not fit the data as well, see Figure 9. However, the Normal distribution 
was the best fit of the distributions tried. Upon studying the patterns of the data with the distribution 
overlay curves, it was obvious that a scaling factor was needed to model some of the data. The model 
with scaling factor that will be used is defined as: 
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1 


e 


( 1 ) 


where 


f(x) = “ 

P 


yf2no 


-(*~t0 2 

2a 2 


a = maximum number of occurrences from the collected voltages 

P = maximum value previously calculated for the Normal distribution 

x = collected voltages 

Ll = mean of collected voltages 

a = standard deviation of collected voltages 


The scaling factor that was used came from multiplying the Normal distribution by the maximum number 
of occurrences from the collected voltages divided by the maximum value calculated for the Normal 
distribution. Figure 10 shows a plot with the original test data, the Normal distribution curve, and the 
Normal distribution curve with a scaling factor incorporated. The curve with the scaling factor fits the 
data and will be used in further analysis. At this point the first goal of the research has been met which 
was to develop and validate a mathematical model of the electromagnetic fields coupling into the 
equipment. 
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Figure 7. Subplot of open-loop test data. 
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Table 5. Distributions Used to Fit the Open-Loop Experiment Data 


Distribution 

Probability density function 

Uniform 

f(x) = — - — , a<x<b; 

b-a 

0, elsewhere 

Normal 

-(x-g) 2 

1 2 

f (x) = — -== — e 2cr , -°°<x<°°; 

V27to 

— °°<x<°°, ct > 0 

Weibull 

f( x ) = -^-x a_1 e _ 

e a 

0, 

f-T 

w , x > 0; a, 0 > 0 
elsewhere 

Rayleigh 

ft \ x Uo 2 
f(x) = — r-e v *■ 

a 

, x > 0 

Beta 

f( x ) = — t-t-t-tx 01 *(l-x)^ 1 , 0<x<l;a, (3 > 0 
v ' r(a)r(p) v ' 

0, elsewhere 

Gamma 

r(r) = j^°x r_1 e _x d x , r>0 

Exponential 

—x 

f(x) = -e * , x>0, 0>O 
0, elsewhere 

Pareto 

Q 

f (x)= 0 < x < °°, 0 < 0 < °° 

(l+x) e+1 
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Figure 8. Open-loop test data with good 
fit of Normal distribution overlays. 
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Figure 9. Open-loop test data with poor 
fit of Normal distribution overlays. 
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Figure 10. Open-loop test data with Normal distribution 
and scaling factor overlays. 


Open-loop test data 
Normal distribution 
Normal distribution with scaling factor 


The focus will now be on the second goal of the research which is to use the model to 
extrapolate/predict measurements that are outside of the given data. The first effort of this statistical 
analysis involved finding the least square error using the equation: 

n 

L = X[ yi -p( Xi )] 2 (2) 

i=l 

where y is the test data and p(x) is the calculated Normal distribution value [16]. 

The method of least squares chooses solutions with coefficients that minimize the sum of the 
squares of the vertical distances from the data points, which are presumed to be polynomial [16]. The 
best-fit polynomial is the one with coefficients that minimize the function L. Figure 1 1 shows a 
representative plot where the least square error was calculated for each field strength and reference 
voltage for one particular frequency and signal. The text in the individual subplots indicates which set of 
parameters best fit the curve. Some data have better least square error estimates than others do, while the 
one that is represented in the subplots is the best fit (lowest value) of the least square error combinations 
tried. 


There were numerous least square error parameters to look at. Table 6 describes the different 
means and variances that were calculated and compared to find the lowest value of the least square error. 
The 'Lxy' number displayed in the subplots represents the combination of parameters used to find the least 
square error. The 'x' (which is 1,2, or3) indicates which reference voltage was used while V (which is 0 
through 15) describes the mean and variance combination used to test different least square errors. 
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Figure 11. Open-loop test data with least square error overlays. 


Table 6. Mean and Variance Combinations Used to Determine the Least Square Error Value 


y 

Mean 

Variance 

0 

mean of matrix 

variance of matrix 

l 

mean of matrix 

sample variance 

2 

mean of matrix 

variance of non-repeated matrix 

3 

midpoint of matrix 

variance of matrix 

4 

midpoint of matrix 

sample variance 

5 

midpoint of matrix 

variance of non-repeated matrix 

6 

midpoint of sorted, non-repeated matrix 

variance of matrix 

7 

midpoint of sorted, non-repeated matrix 

sample variance 

8 

midpoint of sorted, non-repeated matrix 

variance of non-repeated matrix 

9 

mean of non-repeated matrix 

variance of matrix 

10 

mean of non-repeated matrix 

sample variance 

11 

mean of non-repeated matrix 

variance of non-repeated matrix 

12 

voltage with most occurrences 

variance of matrix 

13 

voltage with most occurrences 

sample variance 

14 

voltage with most occurrences 

variance of non-repeated matrix 

15 *amp factor 

maximum voltage 

variance of non-repeated values 


Once the least square error with the lowest value was determined, the mean and variance that 
supported that least square error were saved to different four-dimensional matrices (field strength x 
reference voltage x signal x frequency). Using these new parameters, the values for all of the field 
strengths for one frequency and signal were plotted to graphs. Figure 12 shows a typical graph of the new 
parameters. The top subplot represents the variances collected using the least square error method. Since 
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all of the reference voltages have similar variances, this is a good indicator that the analysis was done 
properly. The bottom subplot represents the means that were collected. This particular signal had 
minimal fluctuation in the means collected. All of the graphs produced were then studied to determine 
what direction to take for further analysis. 
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Figure 12. New mean and variance parameters generated 
from the least square error value. 


The means found using the least square error method varied only slightly for most of the signals, 
so the first focus was on how to extrapolate the variance. This effort involved extrapolating data by hand 
from the new parameters found from the least square error method. This technique produced eight signals 
as having the most variable movement. This hand plotting technique was tedious and became unwieldy 
after a few attempts. The next step was to use an automated computer tool to extrapolate the values to the 
predetermined level of lOOOv/m. The values that were extrapolated depended on the frequency that was 
being examined, as shown in Table 7. 

The basic fitting interface allows the data to be fit using an interpolant or a polynomial (up to 
degree 10). It permitted multiple fits to be plotted simultaneously so that comparison for the given data 
set could be made. It also permitted examination of the numerical residual of a fit, to evaluate the fit, and 
to save the evaluated results to a Matlab variable. The results of the numerous basic fitting interface 
computations were usually cubic in nature, but there were several that were linear or quadratic. 

The best fit was determined by examining the numeric value of the norm of the residuals and the 
residual plots. The fit residuals are defined as the difference between the ordinate data point and the 
resulting fit for each abscissa data point. The norm of the residuals is a measure of the goodness of fit, 
where a smaller value indicates a better fit than a larger value. During this analysis, several of the fits had 
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negative values when extrapolated and it was necessary to return to the basic fitting interface and consider 
a different fit. Once the best fit was determined, the fit was extrapolated using the values from Table 7. 


Table 7. Field Strengths Extrapolated at Each Frequency 


Frequency (MHz) 

Extrapolated Field Strengths (v/m) 

250 

300, 400, 500, 600, 700, 800, 900, 1000 

350 

300, 400, 500, 600, 700, 800, 900, 1000 

450 

650, 700, 750, 800, 850, 900, 950, 1000 

550 

650, 700, 750, 800, 850, 900, 950, 1000 

650 

650, 700, 750, 800, 850, 900, 950, 1000 

750 

650, 700, 750, 800, 850, 900, 950, 1000 

850 

650, 700, 750, 800, 850, 900, 950, 1000 

950 

650, 700, 750, 800, 850, 900, 950, 1000 


Figure 13 shows a typical result from the basic fitting interface tool. The top subplot shows the 
test data, a quadratic fit, and a cubic fit. The two fits are very close and must have another step 
introduced in order to find the best fit. Using the basic fitting interface tool the analysis can be expanded 
to include the residuals, this is represented in the bottom subplot of Figure 13. This subplot shows a plot 
of the residuals with the value of the norm of the residuals. As stated before, the lower the value of the 
norm of the residual, the better the fit. Therefore, in this case, because the norm of the residual for the 
cubic fit was lower than that of the quadratic fit, the cubic fit was used to extrapolate the values of the 
variance. 
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Figure 13. Plot generated by Matlab's basic fitting interface. 
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Figure 14 shows the extrapolated values in the top subplot. The '+' marks show the previously 
calculated variances using the least square error method, while the 'O' marks show the values that were 
extrapolated using the basic fitting interface tool. These values were calculated with the tool and then 
stored. 


hdot@450 MHz 



residuals 



Figure 14. Plot of extrapolated variances using 
Matlab's basic fitting interface. 


Figure 15 shows a plot with the mean and median of the test data plotted. The upper dotted line 
in the figure represents the median while the lower dotted line is the mean. In most cases, the median 
seemed to represent the test data better and was used for the mean value in further analyses. 

Once the variances and means had been obtained, the developed mathematical model could be 
used to plot the Normal distribution using these extrapolated values. Figure 16 shows subplots of the 
scaled Normal distribution using the extrapolated values for the mean and variance. These plots follow 
the same trend that was referred to before, as the field strength increases the data becomes flatter and 
spreads out over a broader voltage range. 

After studying the data and applying numerous statistical tests, including Chi square, student t 
distributions, and confidence intervals, it has been determined that the field strength can be increased up 
to 1000 v/m without damaging the equipment. Therefore, at this point the initial objectives of this 
research have been accomplished. 
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Figure 15. Mean/median computed using 
Matlab's data statistics tool. 
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Figure 16. Normal distribution using the 
extrapolated variance and mean 
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4 Conclusions 


The preliminary objectives of this research have been accomplished. A mathematical model of 
the electromagnetic radiation coupling into the flight control computer has been developed and validated. 

Thousands of representations of the open-loop test data were plotted for comparison. A trend in 
the data was discovered that warranted fitting distribution curves to the data. Because the data had no 
well-defined shape, selected distribution curves were fit to the data and the Normal distribution was found 
to have the best overall fit. It was determined that the Normal distribution needed a scaling factor. Using 
the experiment data, a scaled normal probability density function was developed and validated. Using the 
scaled Normal distribution as the standard, the least square error was then calculated. The mean and 
variance that supported the least square error with the lowest value were plotted for all field strengths for 
one frequency and signal. These plots were used to extrapolate the variances and the means. 

Using the extrapolated values for the variances and means, the developed mathematical model 
was used to plot the Normal distribution for the field strengths that had been extrapolated. These plots 
followed the same trend as the test data: as the field strength increases, the amplitude of the voltage 
flattens and the voltage range covers a broader region. Examining the extrapolated plots led to the 
conclusion that more tests can be performed and the field strength can be increased to the predetermined 
field strength of 1000 v/m without damaging the equipment. 
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Appendix A 

Reference voltages sent for each signal 


Signal 

Min 

Zero 

Max 

cas 

-4.7800 

0 

4.7750 

delac 

-5.0000 

0 

5.0000 

delec 

-5.0000 

0 

5.0000 

delrc 

-5.0000 

0 

5.0000 

deltc 

-5.0000 

0 

5.0000 

epr 

-4.6880 

0 

4.6830 

eta 

-4.3850 

0 

4.3800 

gamma 

-2.8810 

0 

2.8760 

gse 

-2.5000 

0 

2.4980 

hddot 

-2.5000 

0 

2.4980 

hdot 

-5.0000 

0 

4.9950 

phidg 

-5.0000 

0 

4.9950 

phidt 

-4.9950 

0 

4.9900 

psidt 

-4.4970 

0 

4.4920 

qbdg 

-2.5000 

0 

2.4980 

ralt 

-4.9900 

0 

0 

rbdg 

-2.4950 

0 

2.4930 

spl 

-2.9100 

0 

2.9050 

spr 

-2.9100 

0 

2.9050 

tas 

-4.7800 

0 

4.7750 

tbax 

-0.9910 

0 

0.9890 

thrtrm 

-5.0000 

0 

4.9950 

tk 

-4.4970 

0 

4.4920 

vele 

-2.8710 

0 

2.8660 

vein 

-1.2500 

0 

1.2490 

vgs 

-2.5000 

0 

2.4980 

vgsdt 

-5.0000 

0 

4.9950 

ycg 

-1.2490 

0 

1.2480 
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Appendix B 

FORTRAN B737 Autoland Simulator Variables (Partial List) 


variable name 

description 

cas 

calibrated airspeed in knots 

delac 

aileron command in degrees (+ right wing down) 

delec 

elevator command in degrees (+ nose down) 

delrc 

rudder command in degrees (+ yaw left) 

deltc 

throttle command in degrees (always +) 

epr 

engine pressure ratio 

eta 

localizer error in degrees (+ right) 

gamma 

flight path angle in degrees 

gse 

glide slope error in degrees (+ above beam) 

hddot 

vertical acceleration in fps 2 

hdot 

rate-of-climb/descent (+/-) in fps 

phidg 

roll angle in degrees 

phidt 

roll rate, 1-axis in radians/second 

psidt 

yaw rate, 1-axis in radians/second 

qbdg 

pitch rate, body axis in degrees/second 

ralt 

runway altitude in feet 

rbdg 

yaw rate, body axis in degrees/second 

spl 

left spoiler in degrees (+ spoiler up) 

spr 

rt. spoiler in degrees (+ spoiler up) 

tas 

true airspeed in knots 

tbax 

is 1 . if auto stab is driving, else 0. 

thrtrm 

trim throttle position in degrees 

tk 

track angle in degrees 

vele 

east inertial velocity in fps 

vein 

north inertial velocity in fps 

vgs 

ground speed in fps 

vgsdt 

ground acceleration in fps 2 

ycg 

a/c y-position in ft (+ right of center line) 
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